Fidelity optimization for holonomic quantum gates in dissipative environments 
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We analyze the performance of holonomic quantum gates in semiconductor quantum dots, driven 
by ultrafast lasers, under the effect of a dissipative environment. The environment is modeled as a 
thermal bath of oscillators linearly coupled with the electron states of the quantum dot. Standard 
techniques make the problem amenable to a numerical treatment and allow one to determine the 
fidelity as a function of all the relevant physical parameters. As a consequence of our analysis, we 
show that the disturbance of the environment can be (approximately) suppressed and the perfor- 
mance of the gate optimized — provided that the thermal bath is purely superhomic. We conclude 
by showing that such an optimization is impossible for ohmic environments. 

PACS numbers: 03.67.Lx 
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I. INTRODUCTION 

In the last years holonomic quantum computation 
(HQC) has proved to be a viable and promising approach 
to quantum information processing and quantum compu- 
tation P]. According to this approach, quantum informa- 
tion is encoded in an n-fold degenerate eigenspace of a 
family of quantum Hamiltonians depending on dynam- 
ically controllable parameters. Recently, concrete pro- 
posals for quantum computation have been put forward, 
for both Abelian 0] and non- Abelian holonomies 0, 13 • 
The so-constructed gates depends only on global geomet- 
rical feature, e.g., the solid angle spanned in the param- 
eter space, and because of this, it is believed that they 
are robust against errors affecting the physical parame- 
ters controlling the gates themselves (e.g., laser pulses). 
This expectation has been confirmed by recent investiga- 
tions 0. 

An important open problem is whether holonomic 
quantum gates are stable under the effect of the environ- 
ment 0. In this paper we argue that holonomic gates 
have indeed a good performance when the effect of the 
environment is taken into account. We show this on the 
basis of a simple and idealized model which covers situ- 
ations describing electron states in quantum dots rang- 
ing from excitons Q to optically active spin-degenerate 
ones (for a related investigation about environmental 
effects on semiconductor based quantum gates see Ref. 
0). For such gates, the main source of dissipation is 
due to phonons with superohmic spectral density P|. It 
turns out that these gates manifest a rich structure when 
the "control parameters" are changed. By varying in a 
suitable way the adiabatic time, the superohmic effect 
can be minimized and suppressed. For completeness, we 
have extended our analysis to an ohmic environment. 

In Sec. ^ we review the HQC model, introduce a 
model for the dissipative environment and write the mas- 
ter equation to solve. In Sec. IIIII we descibe the com- 
puter simulations and specify the kind of environment 



we consider. From the numerical results we deduce non- 
trivial behavior of the fidelity (the gate performance ex- 
timator) which allows us to minimize the decoherence 
effects of a superohmic environment. We show that anal- 
ogous minimization cannot be done with other kinds of 
environment (ohmic). 



II. HQC IN A DISSIPATIVE ENVIRONMENT 

The physical system we consider is constituted by three 
degenerate (or quasi-degenerate) states (|-f), |— ), and 
|0)) optically connected to another state |G). Every de- 
generate state is separately addressed by polarization or 
frequency selection with a laser. This model describes 
various quantum systems interacting with a laser radi- 
ation field; here, we deal with quantum states in semi- 
conductor quantum dots such as excitons 4] and spin- 
degenerate electron states 7,] • 

The (approximate) Hamiltonian modeling the effect of 
the laser on the system is (for simplicity, /i = 1) 53| 

Ho{t)= E m{j\+e-"-%{t)\j){G\) + h.c., (1) 
j=+-,o 

where flj (t) are the time dependent Rabi frequencies de- 
pending on controllable parameters, such as the phase 
and intensity of the lasers, and e is the energy of the 
degenerate electron states. The Rabi frequencies can 
be modulated within the adiabatic time tad (which co- 
incides with the gating time) to produce a loop in the 
parameter space and thereby realize the periodic condi- 
tion Hoitad) = Ho{0). 

It should be observed that different loops in the pa- 
rameter space produce different holonomic operators. 
Here we consider two different sets of Rabi frequen- 
cies, that is, two different choices of the time depen- 
dent functions flj{t) in Eq. According to the 
holonomic approach, the corresponding unitary evolu- 
tions occurring in the adiabatic time define the unitary 
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FIG. 1: Fidelity for gate 1 as a function of temperature for 
different ks (expressed in (meV)"'^) in the presence of a super- 
ohmic environment. The points are the computer simulation 
results and the curves have the form JT = 1 — X]j=± rijVHad 
ij)- = 3 ■ lO"'^ and r)+ — 0.7). Inset: Fidelity ^ as a func- 
tion of fca for two different temperatures (Ti/Sl = 1.6 ■ 10~^ 
and T2/Q. = 5 ■ 10"^). Parameters: e = 1 eV, f7 = 25 meV, 
tad ~ 7.5 ps and cuc ~ 0.5 meV. 



transformations associated with the holonomic quantum 
gates We shall hereafter refer to the gates defined 
by these two different sets of Rabi frequencies (which 
form a complete set of single qubit gates) as "gate 1" 
and "gate 2"; they correspond, respectively, to the uni- 
tary operators Ui = e^^l+X+l and U2 — e*'5'^« (where 

^. = »(l+)(-|-h>(+l))- 

The Hamiltonian has four eigenstates: two eigen- 

states \B±) with time-dependent eigenvalues e± (called 
bright states) and two eigenstates |£'i,2) with constant 
and degenerate eigenvalue e (called dark states) . To con- 
struct a complete set of holonomic quantum gates it is 
sufficient to restrict Rabi frequencies ilj (t) such that the 
norm of the vectors {n,j{t)} is time independent. Un- 
der this condition, it can be easily shown that the two 
dark states have energy e and the two bright states have 
time independent energies 



e± = (e±y^^T4?^2)/2. 



(2) 



The adiabatic condition is simply ^tad 3* 1- 

The environment is described as an ensemble of har- 
monic oscillators linearly coupled to the system [lQ | , with 
total Hamiltonian 
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H = Ho{t) + J2i 



9 

Pa 

2mn 



1 



+ CaXaA) (3) 



The interaction should break the degeneracy of the 
degenerate states, a condition that is easily fulfilled 
by assuming the operator A to be of the form A = 
diag(0, 1,0, -1) in the basis \G), |+), |0), and |-). (Eq. 
p|) includes implicitly the standard renormalization term 

We now consider the time evolution of the reduced den- 
sity matrix of the system, determined by the Hamiltonian 



FIG. 2: Fidelity for gate 1 subject to superohmic spectral 
density as a function of tad {T/Q = 8.5 ■ 10"^, fcs = 10"^ 
(meV)^^ and Qtad ~ 280). The adiabatic time is normalized 
to 5 ps. Parameters as in Figure^except tad- 



(PJ. We rely on the standard methods of the "master 
equation approach," according to which the effect of the 
environment is considered in the Born approximation, 
and the environment is assumed to be at each time in its 
own thermal equilibrium state at temperature T. One 
has Jj 



dT{g{T)[AA'p{t - r) - A'p{t - r)A] 
+9{-r)[p{t - r)A'A ~ Ap{t - r)A']} (4) 



where, p{t) denotes the time evolution of the reduced 
density matrix of the system in the interaction pic- 
ture, e.g., p{t) = W{t,0)pU{t,0), where U{t,0) = 
T(exp— i J*dt'Ho{t')) and T is the time ordered prod- 
uct. In the above equation A and A' stand for A{t) and 
A{t — t) (again, the superscript tilde means time evolu- 
tion in the interaction picture). In Eq. I^J the effect of 
the environment is included in the function (for simplic- 
ity, fcs = 1) 



9{r) 



J{uj)[coth{-^)cos{LJT) 





i sm{ujT)]dLu (5) 



where the spectral density J(w) is defined in a standard 
way in terms of oscillator parameters in Eq. Q and its 
typical behavior, for physical environments in the low 
frequency regimes, is proportional to w*, with s > 1 j^; 
the asymptotic decay of the real part of ^(r) defines the 
characteristic memory time te of the environment. 

We solve Eq. for an environment with memory 
time Te faster than the time scale of the variation of the 
density matrix td (estimated self-consistcntly) , so that 
p{t — r) « p{t) (Markov approximation). In this approx- 
imation, Eq. Q simplifies and assumes a convenient 
form — for the numerical analysis we wish to perform — in 
the Schroedinger representation: 



p{t)^-^[Ho{t),p{t)]-Cip) 



(6) 
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FIG. 3: Fidelity T for gate 1 subject to superohmic (dashed 
line), ohmic (dot line) and both (solid line) spectral densities 
as a function of T/Q,. Inset: Fidelity T for superohmic plus 
ohmic environment as a function of tad with ^tad ~ 280 (fci — 
4 • 10"'' and ks = 10"^ (meV)"^). Parameters as in Figure Q 



FIG. 4: Fidelity for two qubit gate subject to superohmic and 
ohmic spectral density as a function of T/f2 for fci = 10~* 
and fca = 10~^(meV)~^. Parameters as in Figure except 
tad ~ 0.8 ns and Q — 0.2 meV. 



with 

/•OO 

C{p) - J dT{9{r)[AA"p{t) - A" Pit) A] 

+g{-TMt)A"A^Ap{t)A"]}, (7) 

where A" — W {t ~ t, t)AU {t — t, t) , and, in the adiabatic 
approximation, U{t ~ T,t) w exp(iri/o(i))- As usual, 
the upper extreme of integration is extended to infinity 
because the evolution time is much longer than te- 

III. COMPUTER SIMULATIONS AND 
RESULTS 

In order to estimate how the environment affects the 
performance of the ideal holonomic gates determined 
by the "dissipation-free" Hamiltonian Ho{t), we use the 
standard notion of fidelity 

T = \/ {■llJid{tad)\p(tad)\i'id{tad)) (8) 

where \ipiditad)) is the state in which the initial (pure) 
state evolves, in the adiabatic time tad, under the action 
only of Ho{t), while p{tad) is the solution of Eq. ©, com- 
puted at the same time, and for the same initial (pure) 
state. In order to avoid dependence on the latter we have 
taken a suitable average on the initial states. We make a 
sampling of the initial logical states (combination of the 
logical states |-|-) and |— )) on the Bloch sphere. We add 
the possibility of an error in the preparation of the ini- 
tial state with the population of the non-logical state |0). 
This can be due to the imprecise control of the lasers. 
The initial state has the form a\-\-) -\- /?|— ) -I- ?y|0) with 
\a? + W + \V? = 1 and jryp = 0.1. In the following, 
with a slight abuse of notation, we shall denote by the 
same symbol T this averaged fidelity. 

An essential ingredient of our analysis is the spectral 
density entering in Eq. |SJ|. For the electronic states in 



quantum dots the decoherence effects are principally due 
to phonons. Single phonon processes are described by 
superohmic spectral densities with J(w) = fcacj^e"^"/'^''-* 
The high frequency cut-off ujc is due to the pla- 
nar confinement in the quantum dot. The adimensional 
coupling constant allows the description of different 
kinds phonon-carrier interactions in semiconductor ma- 
terials including deformation potential, piezoelectric and 
spin-orbit p^ . 

The results of numerical simulations for the fidelity J- 
of gate 1, as a functions of the temperature are shown in 
Fig. n] at low temperature the fidelity is constant and 
decreases linearly as the temperature increases. To un- 
derstand how this behavior comes about it is convenient 
to express the dissipative part of C{p) in Eq. © (in the 
dark- bright state basis) as 0: 

^(<5;m ^ '^nrrk + ^nk ^ ^ Irrk ~ ^ Imnk ~ ^tmnk)Pkl ■ 
kl r r 

Here, F^^^. = V^^,Kimnk , with Ki„ink depending only 
on laser parameters and 

^nkM = J(o.„fc)(coth(ig^) T 1) (9) 

are the transition rates between n and k states due to 
phonons. In passing, we note that these are indeed the 
rates that could be guessed by a straightforward applica- 
tion of Fermi's golden rule to the interaction terms of Eq.( 
|5Jl. Finally, observe that the frequencies ujnk = e-n — £k 
represent the energy differences in the dark-bright space, 
i.e., e„ = e, e±, where e± are given by Eq. (0); for fl, 
UJnk = 0,e,f]2/£. 

Equation shows that, with superohmic spectral den- 
sity, the only relevant transition is the one with ujnk — 
ri^/e giving the transition rates r^(f7^/e). In fact, for 
^nk = e the gaussian cut-off with cjc ^ e produces neg- 
ligible rates, and for degenerate dark state (w„fe = 0) 
the rates vanish. These considerations, together with 
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the explicit form of the rates given by Eq. @ provide 
a compelhng explanation of the temperature behavior of 
the fidelity in FigQ] Moreover, we have found that the 
numerical results are fitted by means of the function 

T = \-tadY, ^^r^ (10) 

(where rjj are two real parameters) . This behavior is also 
manifest by considering the fidelity as a function of the 
coupling parameter ^3 (see inset in Fig. 0. A similar 
dependence of the fidelity on the transition rate has been 
found in Ref. 0. Note that in that case the authors 
have only a single transition process (absorption to higher 
states or emission to lower states) while we need to take 
into account both absorption and emission processes for 
the transition to the higher state. This is due to the 
different master equation solved. In fact, they solve a 
strictly second order master equation, while our eq. © is 
"self-consistent" . The similarity in the results is derived 
from the small value of the coupling constant but for 
higher values the two approximations diverge from each 
other and the numerical results cannot be fitted by such 
an elementary function. 

Eq. H10() is particularly important in that it allows 
one to predict how the fidelity behaves when the pa- 
rameters of the system are modified. To this end, 
first of all note that, by keeping the adiabatic pa- 
rameter constant [VLtad — a ^ const), the rates 
become a non-trivial function of the adiabatic time. 
Then, by writing the rates explicitly, one obtains F^ oc 
l/t6^coth(aV(ei2^T))exp(-(aV(t2^e^e))'). Thus, it 
follows from Eq. ^ that the fidelity should have a pro- 
nounced minimum as a function of tad- This behavior 
is confirmed by the computer simulations presented in 
Fig. 13 By varying and tad (e.g., by acting on the 
lasers) the position of the fidelity minimum in Fig. |5]can 
be shifted and the effect of superohmic environment can 
be suppressed. It seems to us that this is an interesting 
result. 

Before discussing the limitations of this result, we 
would like to comment on the approximations upon 
which our analysis relies. As we have already antici- 
pated, in Eq. © the Markov approximation is appro- 
priate when the memory time te is small with respect to 
the time scale of variation of the density matrix td ■ Eq. 
Q leads to the estimate te « l/(27rr); while Eq. 10, 
for a superohmic spectral density and T <C i^^/e, leads 
to T£) — {e/W^)^ /ks. Note that the conditions of validity 
of the Markov approximation, te < td, readily translate 
into a temperature regime, namely T > Tm = k3{Q? / e)^ . 
With our choice of parameters, we have 1.2 • 10"* < 
Tm/^ < 1.2 • 10"'^ (depending on k^ value), which, in 
our simulations, is a very low temperature. 

Since the possibility of suppressing the superohmic ef- 
fects is indeed surprising, one may wonder whether a sim- 
ilar possibility arises for more general environments, e.g. 
for ohmic environments. Though ohmic environments 
are typical of baths of conduction electrons (i| , even for 



phonon baths, which are typically superohmic, it seems 
possible that the spectral density contains an ohmic part; 
this is presumably due to higher order contributions such 
as two phonon processes PJ^. 

Be that as it may, we found it interesting to extend 
our analysis to an environment with the spectral den- 
sity J{uj) = fcitt'e"^'^/'^'^)'' with ki < 1. As is easily 
seen, Eq. @ for ohmic rates leads to completely differ- 
ent results. This is due to the presence, in the ohmic 
case, of transitions between degenerate states which are 
absent in the case of a super homic environment. This 
difference has striking consequences: the transitions be- 
tween degenerate states give contribution to the rates © 
which are linear in T (while the transitions between non- 
degenerate states have the same temperature behavior as 
the superohmic case) . This difference is confirmed by the 
computer simulations in Fig. |31 which shows the fidelity 
for the superohmic and ohmic environments, and the sum 
of the two contributions as a function of T. This curve 
clearly shows that the presence of an ohmic environment 
changes dramatically the fidelity behavior, whence it fol- 
lows the impossibility of extending to the ohmic case the 
results previously obtained changing the adiabatic time. 
This conclusion becomes very clear if one compares Fig. 
12 with the inset in Fig. |3| It is not possible anymore to 
optimize the fidelity by changing the parameters. 

Before concluding, we briefly mention three points (see 
[1^ for a thorough discussion). First, the computer sim- 
ulation for gate 2 confirms the results found for gate 1. 
Second, our analysis extends (almost straightforwardly) 
to the two-qubit gate proposed in Ref. 4J: Fig. ^ 
shows the preliminary results for the two-qubit gate (for 
a non random initial state) as a function of temperature 
and coupling constant for superohmic and ohmic environ- 
ments; the behavior of the fidelity is analogous to that of 
the one qubit case (Fig.|31). Third, for the ohmic case, a 
careful study of the (relatively) high temperature behav- 
ior shows the limitations of the Lindblad approximation 
for the reduced dynamics. 



IV. CONCLUSIONS 

To sum up, the upshot of our analysis is twofold. The 
good news is that it is possible to optimize the fidelity 
for the kind of environment which is usually considered 
for electron states in quantum dots — a superohmic en- 
vironment caused by electron-phonon interactions. The 
bad news is that such optimization does not go through 
an ohmic environment, e.g., produced by the same su- 
perohmic phonon bath through two phonon processes. 
Thus, particular attention should be paid in modeling 
the environment, since the presence of a weak ohmic en- 
vironment dramatically changes the holonomic gate per- 
formance. For these reasons it is crucial to dispose of 
experimental investigations on the nature of the envi- 
ronmental spectral densities in semiconductor quantum 
dots. 
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